Probing isotope effects in chemical reactions using single ions 
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Isotope effects in reactions between Mg^ in the 3p^P3/2 excited state and molecular hydrogen 
at thermal energies are studied through single reaction events. From only ~250 reactions with 
HD, the branching ratio between formation of MgD^ and MgH^ is found to be larger than 5. From 
additional 65 reactions with H2 and D2 we find that the overall decay probability of the intermediate 
MgHj, MgHD^ or MgDj complexes is the same. Our study shows that few single ion reactions 
can provide quantitative information on ion-neutral reactions. Hence, the method is well-suited for 
reaction studies involving rare species, e.g., rare isotopes or short-lived unstable elements. 
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Isotope effects often play an important role for the out- 
come of chemical reactions. For instance the chemical 
composition of interstellar clouds is strongly influenced 
by isotope effects in certain reactions [1]. In laboratory 
experiments, isotope effects observed in isotopic analogs 
of chemical reactions can provide details about the reac- 
tion dynamics. Substitution reactions of the type F+HD 
and Cl+HD are among the simplest reactions where iso- 
tope effects can be present. They were already studied 
when the first chemical lasers were developed, in order 
to understand the population inversion mechanism and 
to identify the laser transitions [2]. These studies to- 
gether with studies of the F+D2 reaction strongly stimu- 
lated the whole field of reactive scattering [3] . Especially 
the resonance effects observed in the F+H2 reaction and 
isotopic analogs [4] have been subject to numerous ex- 
perimental and theoretical studies finally resulting in a 
much iniproved understanding of this benchmark reac- 
tion [1, m. In another series of experiments, reactions 
between a beam of ground state atomic ions and H2, HD 
and D2 have been studied. Strong isotope effects have 
been observed and, e.g., in reactions between alkaline 
earth ions (X+) and HD, it was found that XD+ forma- 
tion is preferred for some alkaline earth ions, while XH+ 
formation is preferred for others [6, 7]. 

Within the last few years, new techniques have 
emerged for ion- neutral reaction studies, e.g., of cold re- 
active collisions (sj , single ion reactions [§] as well as reac- 
tion dynamics by applying crossed molecular beam imag- 
ing [lO|. These techniques open up new possibilities in 
the field of ion-neutral reactive scattering. 

In this Letter, we investigate isotope effects in reac- 
tions of single ions with the isotopologues of molecu- 
lar hydrogen, a model system for ion-neutral reactions. 
More specifically, we consider reactions at thermal ener- 
gies with Mg^ in the 3p^P3/2 excited state (excitation 
energy of 4.4 eV). Due to the simple internal structure of 
the reaction partners, the studied reactions represent a 
simple test case for reaction dynamics involving an elec- 
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FIG. 1: Sketch of the experimental setup for non-destructive 
identification of molecular ions. See text for details. 



tronically excited atomic collision partner. We make use 
of an experimental technique with almost 100% efficiency 
in detecting single reaction events [9]. With a total of 
only about 300 reactions, the branching ratio between 
the reactions 

Mg+(3p 2P3/2) + HD ^ MgD+ + H (1) 
^ MgH+ + D (2) 

has been found to be larger than 5. In reactions with 
H2, HD and D2 the decay paths leading to either MgH^ 
or MgD^ formation from a MgH2", MgHD+ or MgD2" 
complex have been found to be equally likely within sta- 
tistical uncertainties. Our experiments demonstrate the 
prospects for similar single molecular ion studies, e.g., 
using state prepared molecular ions [11], more complex 
molecular ions [12I , or of astrophysically relevant re- 
actions jl3, EH. The high detection efficiency can fur- 
thermore be useful for studies of reactions involving ions 
of rare species, e.g., superheavy elements [16j . 

In our study we use a linear Paul trap setup which is 
described in detail in Ref. (13. Briefly, as shown in Fig.[Tl 
the trap consists of four cylindrical rods, each sectioned 
into three parts. By applying suitable ac and dc volt- 
ages (not indicated in Fig.[T]) a harmonic confining poten- 
tial is created with oscillation frequencies ujx ^ > 



along the x-, y- and z-axes, respectively. ^^Mg^ ions are 
loaded into the trap by crossing an effusive beam of Mg 
atoms with a laser beam at 285 nm in the trap center for 
resonance-enhanced isotope selective two-photon ioniza- 
tion of ^^Mg The ^^Mg+ ions are Doppler laser 
cooled on the 3s ^S'i/2-3p ^P3/2 transition near 280 nm. 
Individual ions are observed by imaging light, emitted 
spontaneously during the laser cooling process, onto a 
charge-coupled device (CCD) camera. 

Reactions with HD, H2 or D2 molecules are investi- 
gated by first leaking the gasses into the trap chamber un- 
til a steady-state pressure of about 10~^ Torr is reached, 
and then loading two ^^Mg^ ions into the trap. Reac- 
tions exclusively take place with ^^Mg^ ions excited to 
the 3p^P3/2 state in the laser cooling process, since re- 
actions between ^^Mg^ ions in the 3s ground state and 
thermally excited hydrogen molecules are not energeti- 
cally allowed [lO]. After a reaction, the formed molec- 
ular ion stays trapped since its acquired kinetic energy 
only amounts to a small fraction of the trap depth (~1 
eV). Within tens of milliseconds after a single molecu- 
lar ion is formed, it is sympathetically cooled through 
the Coulomb interaction with the remaining laser cooled 
26]\/[g+ form a two- ion Coulomb crystal. Although 

the molecular ion does not emit light, its presence is evi- 
dent from the fact that the remaining laser cooled ^^Mg^ 
ion is located at one of the two initial ^^Mg^ ion positions 
(see Fig.Ej). 

For unambiguous identification of the molecular ion 
species, its mass is determined by applying the identifi- 
cation technique demonstrated for CaO+ ions in Ref. (9|. 
In short, the method relies on exciting the center-of-mass 
(CM) mode of the cold two-ion system along the z-axis 
by applying a voltage oscillating at a frequency Vrnod to 
two of the end-electrodes as shown in Fig. [H When Vrnod 
is equal to the eigenfrequency of the CM mode, the mo- 
tion of the ions is resonantly excited. This excitation is 
clearly visible in the CCD-images as a smearing-out of 
the fluorescence light from the ^^Mg^ ion along the z- 
axis (see Fig. [2]) due to the long exposure time (100 ms) 
compared to the CM mode oscillation period (typically 
~ 10 /is). Since the CM eigenfrequency depends on the 
masses of both ions, the mass of the reaction product ion 
can now be determined [9]. 

In practice, the modulation frequency is scanned in five 
steps through a 100 Hz narrow frequency interval around 
z^rn, where Um denotes the CM mode eigenfrequency for 
one ^^Mg^ ion and one singly-charged ion of mass m, 
and at each step a CCD image is recorded. Such images 
are shown in Fig. [2] for one ^^Mg^ ion trapped simul- 
taneously with another ^^Mg^ ion, a ^^MgH^ ion and 
a ^^MgD^ ion while modulation frequencies of U26^ ^27 
and 1^28 are applied. In the experiments, modulation fre- 
quencies z^28, ^27, ^26, ^25 and 1^24 are applied repeatedly 
until one ^^Mg^ ion has reacted and the reaction prod- 
uct has been identified. This product ion was ^^MgD^ 
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FIG. 2: CCD images of two- ion Coulomb crystals with modu- 
lation applied at frequencies z/26, ^27 and z/28 for identification 
of ^^MgH+ and 2^MgD+ ions. The dashed vertical lines indi- 
cate the position of the trap center along the z-axis. 



or ^^MgH^ in all but two cases where a ^^Mg+ ion pro- 
duced in a charge exchange collision was observed (as in 
Ref. [19] for calcium). The choice of working with the less 
abundant ^^Mg^ ion instead of the dominant ^^Mg iso- 
tope (80% natural abundance) was made to avoid prob- 
lems in distinguishing ^^MgH+ and ^^MgD+ ions from 
25Mg+ and 26Mg+ ions, respectively, formed in charge 
exchange collisions with background gas Mg atoms [19]. 
At the applied molecular gas pressure (~ 10~^ Torr), the 
Mg^ -molecule reaction rate is less than one per minute, 
which means that the above identification procedure can 
in general be applied before both of the initially loaded 
atomic ions have reacted. After a reaction product has 
been identified the trap is emptied, two new ^^Mg^ ions 
are loaded and the experiment repeated. 

In our study of the branching ratio between reactions 
^ and ([2]), the background pressure of H2 is accounted 
for by measuring the number of ^^MgD^ ions, A^MgD+? 
and ^^MgH^ ions, formed at four different ratios 

of the partial pressure of HD and Il2 , Phd / P1I2 • The ob- 
tained results are presented in Fig.[3]and show an increase 
of MgD^ formation with increased HD pressure. 

To understand quantitatively the branching ratio be- 
tween reactions ([T]) and ([2j), we model the reactions as 
two-step processes. In the first step the neutral molecule 
is assumed to be captured by the ^^Mg^ ion at long range 
to form a MgH^ or MgHD+ collision complex and, in 
the second step, a stable MgH^ or MgD^ ion is formed 
from this complex at short range. The capture cross- 
section for the first step is for the present collision en- 
ergies well approximated by the Langevin capture cross 
section a = {e/2eo)y^a/ {fiv'^) [21], where a is the polar- 
izability of the neutral molecule, /i is the reduced mass, 
and V is the relative velocity in the center-of-mass system. 
Since the polarizability of H2 and HD is equal within 
1% [22], in our model we assume a to be identical for H2 
and HD. Due to the low ion oscillation frequencies in the 
trap, the reduced masses are not significantly different 
from the free-free reaction cases. 

The two-step model is valid since the capture range in 
the first step is of the order of (ct/tt)^/^ ~ 6 A which is 
much larger than the extension of the hydrogen molecule 
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FIG. 3: The ratio between the number of formed ^^MgD^ 
and ^^MgH+ ions vs. the relative pressure of HD and H2 gas. 
From left to right the data points correspond to A^MgD+ =17, 
21, 33, and 32 and A^MgH+ =75, 45, 14, and 10. The er- 
ror bars represent statistical uncertainties. A fit of Eq. (|3|) 
to the data results in ^MgD+/''7MgH+ ~ ^'"^^ ^ ^•-'■'^ 
^MgH+/^MgH+ = O-O^lo.Je (black solid line). The red dashed 
line represents Eq. (|3]) in the absence of isotope effects, i.e. 

with ^MgD+/^MgH + 
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and the Mg^ ion as well as the de Broglie wavelength 
relevant for the reaction. If the probability for molecule 
formation in the second step is denoted by 7^, the molecu- 
lar ion formation rate F is then given by Pexc xnxvxcrxr]^ 
where Pexc is the ^^Mg^ excitation probability and n is 
the neutral molecule density. Using that density is pro- 
portional to pressure we find that 
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where 7?^; is the efficiency of M+ (MgH+ or MgD+) 
formation after reaction with M' (H2 or HD). This ex- 
pression is independent of Pexc and of the absolute par- 
tial pressures, which greatly reduces systematic errors. 
The relative partial pressures are measured with a rest- 
gas analyzer and for all values of the ratio Phd/^H2 
we estimate its systematic uncertainty to be less than 
25%. From the parameters determined from the fit (see 
Fig. [3]), the branching ratio ^M^D+/^MgH+ 
ter mined. Taking the statisticaf uncertainties into ac- 
count, we find that the branching ratio is unbound from 
above and has a lower limit of 5, thus demonstrating a 
dramatic intramolecular isotope effect. In addition, the 
ratio (/?MgH+ + ^MgD+)/^MgH+ consistent with unity 
which indicates that the probability of forming a mag- 
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Reaction coordinate 

FIG. 4: Sketch of relevant potential surfaces for the Mg^+HD 
reaction [i^, [IS- On the left (right) hand side the potential 
surfaces are labelled according to C2V (Cs) symmetry. As dis- 
cussed in the text, the reaction probably proceeds by insertion 
of Mg^ into the HD bond on the 1^B2 potential surface, fol- 
lowed by Mg^-D or Mg^-H bond formation. 



nesium hydride ion is equally large from any of the two 
complexes MgHj or MgHD+. 

In a second series of measurements we studied reac- 
tions between Mg^ and a mixture of 43±2% H2 and 
57±2% D2 gas. We observed the production of 40 MgH+ 
ions and 25 MgD^ ions. Using the same two-step model 
as above we find that ^MgH+/^MgD+ =1-5^ O-^- 

Substitution reactions between a ^^MgH^ or ^^MgD^ 
ion and a H2, HD or D2 molecule could potentially give 
rise to a systematic error. Such reactions were, how- 
ever, only observed on two occasions (one MgH^ + D2 
MgD+ + H + D, one MgD+ + H2 ^ MgH+ + H + D) 
in a period of more than 20 minutes and hence they do 
not give rise to systematic errors on the results presented 
here. 

The most striking of the above results is the strong 
intramolecular isotope effect in reactions ([1]) and ([2]). 
This finding cannot be explained by a simple statisti- 
cal model based on an assumption of an equal prob- 
ability for populating energetically accessible states of 
^^MgH^ and ^^MgD^. This assumption only gives rise 
to / ^MgH+ ^ ^ therefore we attribute the ob- 

servecf isotope effect to a dynamical mechanism. In the 
ion beam experiments of Ref. a similar isotope effect 
has been observed in reactions between ground state Mg^ 
ions and HD molecules at center-of-mass energies up to 
11 eV. This was rationalized in terms of an impulsive 
interaction with a thermodynamic threshold. 

A schematic view of the potential surfaces involved in 
the reaction is shown in Fig. [H The analogous potential 
surfaces for MgD2" was explored in a photofragmenta- 
tion study where MgD^ formation was observed after 
laser excitation from the 1 '^Ai state to the red of the 
Mg+(3p)+D2 asymptote From the observed MgD+ 
spectrum it was argued that MgD^ is formed by direct 
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and fast reactions on the 1 surface in C2v geome- 
try through a bond-stretch mechanism as weh as from 
the 1 ^Bi state, possibly through a couphng to the 1 
state. On the 1 potential surface the Mg^ ion be- 
comes inserted in the D2 bond such that the D-D bond 
is stretched and eventually broken and a Mg^-D bond 
is formed [l^, |25|. In addition, in this study no dissoci- 
ation into Mg^(3s)+D2 was observed. In our two-step 
model this corresponds to the values of ^^gH+' ^MgD+ 
and (^^gD+ ^ ^MgH+) being close to unity, which indeed 
is in good agreement with the Mg^+H2 reaction rate 
measured in a previous study 13, Q . 

The investigation of photofragmentation indicate that 
the Mg^+HD reaction discussed in the present pa- 
per proceeds via the 1 ^^2 surface through a bond- 
stretch mechanism that eventually favors the formation 
of MgD+ [2I, [25|. To fully understand the transition 
from a MgHD+ complex to a potential surface favor- 
ing the MgD^+H asymptote rather than the MgH^+D 
asymptote requires a detailed theoretical study. It might 
be necessary to consider the details of the conical inter- 
section which arises from the crossing of the l^Ai and 
1^B2 potential surfaces. Non-adiabatic couplings at the 
conical intersection could give rise to a preference of the 
MgD^ channel over the MgH^ channel. The same mech- 
anism could be responsible for the isotope effect observed 
in reactions with ground state Mg^ ions [6]. 

In conclusion, we have found that reactions between 
Mg^ in the 3p^P3/2 excited state and HD molecules 
at thermal energies preferentially leads to the forma- 
tion of MgD^ rather than MgH^ with a branching ratio 
larger than 5. Additional reactions with H2 and D2 have 
shown that after a reaction complex has formed through 
a Langevin capture process, the molecular ion formation 
efficiencies after capture of a HD, H2 or D2 molecule are 
equal within statistical uncertainties. The efficiencies are 
consistent with unity which is in agreement with obser- 
vations in a previous study of MgD2" photofragmenta- 
tion [25^. Our measurements demonstrate that it is pos- 
sible to determine quantitatively the branching ratios and 
relative reaction rates in ion-neutral reactions by observ- 
ing only a few hundred single reactions. The method 
should in the future be applicable to a variety of studies, 
e.g., of astrophysically relevant reactions, chemistry of 
superheavy elements and reactions involving more com- 
plex molecules. 
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